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—i- ' The equilibrium superfluid phase of 3 He impregnated into high-porosity sil- 

[~-». , ica aerogels appears to be a non- equal- spin-pairing state in zero field at all 

^^ ' pressures, which is generally assumed to be the Balian-Werthamer (BW) 

phase modified by the depairing effects of the aerogel structure. The nu- 
clear magnetic susceptibility played a key role in identifying the B-phase of 
pure 3 He with the BW state. We report theoretical calculations of the nu- 
clear magnetic susceptibility for the BW model of superfluid 3 He in aerogel 
within the framework of the Fermi-liquid theory of superfluid 3 He. Scat- 
*^ ' tering of quasiparticles by the aerogel, in addition to Fermi-liquid exchange 

corrections, leads to substantial changes in the susceptibility of the BW 
phase. The increase in the magnetic susceptibility of 3 He-aerogel compared 
to pure 3 He-B is related to the polarizability of the gapless excitations and 
the impurity-induced local field. The limited data that is available is in 
^A , rough agreement with theoretical predictions. Future measurements could 

;_, ' prove important for a more definitive identification of the ordered phase, 

as well as for refining the theoretical model for the effects of disorder and 
scattering on the properties of superfluid 3 He. 

1. Introduction 

Early torsional oscillator and NMR experiments on 3 He impregnated into silica 
aerogel provided the first evidence for a superfluid transition for 3 He in a disor- 
dered medium. The transition temperature is suppressed in 98% porous aerogel to 
T c /T^ ulk ~ O.7.B0 Sizeable suppressions of the superfluid density and transverse 
NMR shift were also reported. NMR measurements of the susceptibility indicated 
that an equal-spin-pairing (ESP) state was stable at pressures as low as 12 bar,cl 
suggesting that the equilibrium phase of 3 He in aerogel at low pressures and non- 
zero field was not the Balian-Werthamer (BW) state that characterizes the bulk 
B-phase of superfluid 3 He. However, with the addition of roughly four monolayers 
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coverage of 4 He on the aerogel strands, the magnetization was observed tojiecrease 
below the aerogel transition indicative of a non-ESP, B-like, pairing stateE 

More recent measurements have clarified the phase diagram of 3 He in aero- 
gel. Torsional oscillator measurements over the full pressure range (0 — 34 bar) for 
temperatures as low as 0.1 mK provide convincing evidence of a zero-temperature, 
normal-superfluid .phase transition as a function of pressure with a critical pres- 
sure of P c ~ 6 bar. a Additional evidence for a B-like phase over the pressure range 
P = 1.5 — 29.3 bar was reported by Alles, et alJa based on the textural lineshapa 
analysis of NMR data on 3 He in 98% and 99% porous aerogels. Barker et al.cl 
also reported the suppression of the magnetic susceptibility with two monolayers of 
4 He added to suppress the Curie component of the magnetization of solid 3 He that 
plates to the aerogel strands. Their NMR measurements provide evidence for an 
A-like phase very close to T c and a B-like (non-ESP) phase at slightly lower temper- 
atures, with a first-order 'AB' transition exhibiting substantial supercooling. The 
AB transition in 3 He-aerogel was also observed at low, pressures (P == 4.8 bar) in 
relatively high fields (H AB ~ 0.2 T) by Brussaard, et al.Q Gervais, et al.El have shown 
that sound attenuation measurements at frequencies of ~ 15 MHz provide a pow- 
erful tool to map the phase diagram of 3 He-aerogel as a function of magnetic field, 
temperature and pressure. Their recent measurements of T c and T AB as a function 
of field provide strong evidence for a stable B-like phase in the zero-field limit for 
pressures of 25 and 33 bar .Em The AB transition is suppressed below the zero-field 
T c quadratically with field: (T AB {H) - T c )/T c ~ -0.05iJ 2 kG -2 at P = 25 bar. 
These authors also find that the A-like phase exhibits strong supercooling even in 
low magnetic fields, which is the likely explanation for earlier observations of an 
ESP pairing state in the susceptibility measurements by Sprague et al.El 
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Fig. 1. Sketch of the phase diagram of superfmid 3 He in 98% porous aerogel. 
The B-like phase extends to the highest pressures. The quantum phase tran- 
sition is shown at P c — 6 bar, and the field-induced AB transition, Tab{H), 
is indicated by the dash-dotted line. Supercooled A-like 3 He extends down 
to T AB . The bulk transition lines are also indicated. 
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The resulting phase diagram of superfluid 3 He-aerogel (see Fig. Q) differs signif- 
icantly from that of bulk superfluid 3 He. There is a zero-temperature, "quantum" 
phase transition at P c ~ 6 bar separating a disordered normal Fermi-liquid phase 
from a B-like superfluid phase. The A-like phase appears to be an equilibrium 
phase only in the presence of a magnetic held. If there is a region of stable A- 
phase in zero-held then it is. confined to a narrow temperature window below T c 
with |T AB (0) — T c \ < 20/iK.cl The other significant feature of the AB transition in 
3 He-aerogel is the degree to which the A-like phase can be supercooled for pressures 
P-20 bar. 

The non-ESP 'B-like phase' in aerogel is generally assumed to be the Balian- 
Werthamer phase, modified by the aerogel structure. This identification is consis- 
tent with theoretical expectations based on free energy calculations for the possible 
phases of superfluid 3 He-aerogel within the homogeneous scattering modcl.EJ Ho- 
mogeneous, isotropic scattering leads to increased stabilization of the BW phase 
relative to the Anderson-Brinkman-Morel (ABM) phase, the planar phase and the 
polar phase. Nevertheless, the identification of the equilibrium phases of superfluid 
3 He in aerogel has not been rigorously established as it has been in bulk 3 He. Iden- 
tification of the equilibrium phases of 3 He-aerogel is complicated by uncertainties 
in the homogeneity and variations in the microstructure of aerogels of the same 
density, the inhomogeneities induced into the order parameter by pairbreaking ef- 
fects of the aerogel structure, and the increased difficulty in carrying out accurate 
theoretical calculations of the properties of 3 He-aerogel. 



2. 3 He-Aerogel Scattering Model 

In the homogeneous scattering model (HSMLthe aerogel medium is represented 
by a random distribution of scattering centersJiiJ For 98% porosity the typical di- 
ameter of the aerogel strands is d ~ 30 A and the mean distance between strands 
is D ~ 320 A, which is very large compared to the Fermi wavelength, kj 1 ~ A, but 
comparable to or less than the bulk coherence length, £o = fovf /2-kT c q ~ 200 — 800 A 
over the pressure range P = — 34 bar. Thus, the aerogel does not significantly 
modify the bulk properties of normal 3 He, beyond the formation of a couple of 
atomic layers of solid-like 3 He on the silica strands. The dominant effect of the 
aerogel structure is to scatter 3 He quasiparticles moving at the Fermi velocity. 
Such scattering has dramatic effects on the formation and properties of the su- 
perfluid phases. If the coherence length (pair size) is sufficiently long compared to 
the typical distance between scattering centers, then the aerogel is well described 
by a homogeneous, isotropic scattering medium with a mean-free path determined 
by the aerogel geometry. 

Figure |] shows the dependence of the Cooper pair size, £ = fi.w//27rT c , on 
pressure in superfluid 3 He-aerogel calculated for aerogel with a mean free path of 
£ = 1700 A. The HSM may be expected to provide a reasonable first approxi- 
mation to the properties of superfluid 3 He-aerogel, particularly at pressures below 
15 — 20 bar. Indeed the HSM accounts semi-quantitatively for the reduction of T c , 
including the quantum ^ritical pressure P c , and the pair-breaking suppression of 
the order parametcr.EiEj However, the HSM is expected to become a poorer de- 
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Fig. 2. The coherence length of superfluid 3 He in 98% porous aerogel as 
a function of pressure calculated within the HSM for a mean-free path of 
I — 1700 A. D ~ 320 A is the typical distance between scattering centers. 



scription of 3 He-aerogel at higher porosities and higher pressures where the pair 
size becomes comparable to, or smaller than, the mean distance between scattering 
centers. This breakdown of the HSM is evident in the quantitative discrepancies in. 
the pressure dependence of T c and p s /p, particularly for higher porosity aerogels.Lil 
More elaborate scattering models which take into account the similar scales for the 
pair correlation length and the aerogel correlation length, provide more accurate 
scaling behavior for the physical properties of the superfluid phases with £o(-P)Jij 
In the following we calculate the pairbreaking effects of quasiparticle scattering 
by the aerogel, within the HSM, on the nuclear magnetic susceptibility of the Balian- 
Werthamer phase of superfluid 3 He in aerogel. Comparison with existing data for 
the magnetization is also presented. Extensions of this theory to include aerogel 
correlation effects and inhomogeneities in the aerogel medium can also be carried 
out if future measurements warrant. 

3. Magnetization of Superfluid 3 He-B-aerogel 

The nuclear magnetization, M, of normal liquid 3 He at temperatures, ksT <C 
Ef, and fields, jhH <C Ef, is given in terms of the (single-spin) quasiparticle density 
of states at the Fermi level, Nf, the nuclear gyromagnetic ratio for 3 He, 7, and the 
exchange enhancement of the local field given in terms of the Landau interaction 
parameter, Ftf, 

XN=M/H=^, (1) 

where /i = 7/1/2 is the nuclear magnetic moment of the 3 He nucleus; xn is the 
nuclear spin susceptibility of the normal Fermi liquid. 

The effect of the aerogel on the magnetization of the normal liquid phase of 
3 He is expected to negligible. However, the aerogel structure is known to be cov- 
ered with one or two layers of localized 3 He atoms. These surface solid layers 
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contribute a Curie-like susceptibility that obscures the Fermi-liquid contribution 
at low temperatures.0 The surface contribution can be suppressed by the addition 
of a couple of monolayers of 4 He which preferentially plates the aerogel structure. 
The net effect is two-fold: (1) the surface Curie susceptibility is suppressed and (2) 
spin-spin scattering between 3 He quasiparticles and the surface spins is suppressed. 
The cross-section of the aerogel may also be modified by 4 He preplating, but we 
expect this effect to be relatively small. Measurements of the B-like suppression 
of the susceptibility in 3 He-aerogel have so far been reported -Only for two or more 
monolayers of 4 He added to displace the solid layer of 3 He.B'H In the following we 
assume 4 He coats the aerogel surface and consider only nonmagnetic scattering of 
3 He quasiparticles off the aerogel structure. 

The nuclear spin susceptibility of pure superfluid 3 He-B .agrees quantitatively to 
leading order in T c /Ef with the result of Serene and Raineili3 for the susceptibility 
of the Balian-Werthamer state, 

"2 +r( l + l^ 



(l + ^o a ) L 

XB/XN i + *?(§ + Jm^($ + §m^or ' (2) 

where Y(T) is the well-known Yosida function, 

F(T) = 1-ttTV ^—^, (3) 

and A(T) is the B-phase gap amplitude. Indeed this result, and its generalization 
to include nonlinear field .-corrections arising from the pair-breaking effect of the 
nuclear Zeeman energy Jlil'tLa has been used to obtain the £ = 2 Landau parameter, 
F% , from measuxements of the susceptibility and gap distortion of the B-phasc 
collective modes£3 

Below the superfluid transition in aerogel the magnetization is given by 

M = ^( H+ ^)< (4) 

where m is related to a Fermi surface average of the change in spin polarization of 
quasiparticles at the point p on the Fermi surface, 



*rW 



^g(p;e«)- (5) 



The change in spin polarization leads to a reduction in the magnetization; 
m/2^ = — D(T)H, where D(T) is the net depolarization below T c . The change in 
the spin-polarization is determined by competing effects of pairing correlations of 
quasiparticles with S z = 0, and pairbreaking induced by scattering off the aerogel 
structure. The effects of depairing of the S z — Cooper pairs by scattering from 
the aerogel are expected to be significant and to lead to sizeable changes in the 
magnetization. An obvious effect of scattering by the aerogel on the magnetization 
compared to that of bulk 3 He-B is the suppression of T c relative to T c q. However, 
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the magnitude of the susceptibility, particularly at low temperatures, is sensitive to 
the density of quasiparticle states below the gap, e < A, produced by scattering off 
the aerogel. 

In the HSM for isotropic scattering in the Born and unitarity scattering limits 
the generic form of Eq. (0) for the B-phase susceptibility is preserved with the 
replacement of the gap and Yosida functions by impurity-renormalized gap and 
response functions. The results can be summarized by Eq. (0) with the replacement 
of Y(T) — > Y(T) in the Born and unitarity scattering limits, 
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A 2 1 
Y = 1-ttTV — -<^ = } (unitarity), (7) 

where Tjv is related to the mean- free path (£), or scattering rate (1/t), for normal- 
state 3 He quasiparticles by 

r - h - hVf (r\ 

Tn -2^-^1' (8) 

and the gap and renormalized Matsubara frequencies are defined in Eqs. (Eoh and 
(pq). The terms that depend explicitly on the scattering rate arise from impurity 
vertex corrections to the local field. The impurity field, which has opposite signs in 
the Born and Unitary limits, provides a small, but non-negligible, correction to the 
depolarization at low temperatures. These results are derived in the next section, 
and discussed quantitatively in Sec. (||). 

4. Fermi-Liquid Theory for Superfluid 3 He-aerogel 

The properties of superfluid 3 He in aerogel cau.be calculated within the frame- 
work of the Fermi-liquid theory of superiluid 3 He.E 3 J The effects of scattering by the 
aerogel are included within the HSM.Ej The central equation of the Fermi-liquid 
theory of superfluid 3 He is the quasiclassical transport equation, 

ie n r 3 - A-S imp -S FL -v cxt , g(p,R;e„) + ivf ■ Vg = , (9) 

where e n are the Fermion Matsubara frequencies, v/ = w/p is the Fermi velocity 
for excitations near the point p on the Fermi surface and g is the quasiclassical 
Green's function, which is the 4x4 Nambu matrix in particle-hole/spin space. For 
a full discussion of the quasiclassical theory of superfluidJiHe, including the matrix 
notation used here, see the review by Serene and Rainer.tj The matrix structure of 
the propagator in particle-hole space is given by 

g + g ■ <r icr 2 f + icr(T2 ■ f \ , .. 

ier 2 / + f- i<J2<r g + io 2 g ■ <x ia 2 J 

where <Xj (i = 1,2,3) are the Pauli spin matrices, g and g are the spin scalar and 
vector components of the diagonal quasiclassical propagator, and / and f are the off- 
diagonal spin-singlet and triplet pair amplitudes. The diagonal Matsubara Green's 
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function determines the density of states, magnetization, local exchange field etc. 
Fundamental symmetry relations connect the particle and hole components of the 
propagators and self-energies, 

£<J>, £„) = g{-p, -£„) , g(p, £„) = g(-p, -£„) , (11) 

/(p,£„) =/(-p,£„)* , f (p, £„) = -f (-p, £„)* . (12) 

The gradient term in Eq. (W determines the range for spatial variations of 
the pair amplitude and quasiparticle excitations in the superfiuid phases. In the 
following we consider homogeneous equilibrium states, in which case we drop the 
gradient term in the transport equation. The matrix terms in Eq. (M) represent 
different physical contributions to the self-energy and coupling to external fields. 
The self-energy terms can be classified in terms of the expansion—parameters of 
Fermi liquid theory, e.g. k B T c /Ef, l/fc/£ , h/rEf, ^H/Ej, etc.EJ The leading 
order contributions represent the off-diagonal pairing self energy, A, the self-energy 
resulting from multiple scattering of quasiparticles by the aerogel, S imp , the Landau 
molecular field energies, S FL , and the Zeeman coupling of the nuclear moments to 
an external magnetic field, w oxt . 

The off-diagonal self energy defines the order parameter; in matrix form, 

A(P) -1 V A ) ~ \w y <r > £(p) ) ' [16 > 

where we consider pure p-wave, spin-triplet pairing, which is sufficient for investi- 
gating the linear response to the nuclear Zeeman field.Ej The order parameter A 
satisfies the weak-coupling gap equation, 

|e„l<e 



v- / drip t 



A(p) = T ^ J _^( p .p')f(p';£„), (14) 

where the pairing interaction in the spin-triplet, p-wave channel is given by, V f = 
3Vip-p'. The interaction, V\, and cutoff, e c , that enter Eq. ( |14| ) can be eliminated in 
favor of the measured bulk transition temperature using the linearized equilibrium 
gap equation for pure 3 He in zero field (see below). 

In the HSM the scattering of 3 He quasiparticles off the aerogel structure is 
modeled by a random distribution of scattering centers ( "impurities" ) . The impu- 
rity self-energy to leading order in H/tEj is determined by the t-matrix for multiple 
scattering by a single impurity and the mean density of impurities, 

S ,mp (p; £„) = n imp £(p, p; £„) , (15) 

where n imp is the impurity density and t is obtained from the self-consistent solution 
to the t-matrix equation, 

-^- u(t>, P") S(P"; en) i(f>", P'; En) , (16) 

the quasiclassical transport equation for g and the constitutive equations for the 
self energies. The general structure of the impurity self-energy matrix is given by 

a ,- e \ _ ( S lmp + h lmp ■ a A imp io y + A lmp ■ icray \ ( r 

S toP (P, £n) - (^ iay + Ai ^ . l(jy(T s ^ + ^^ . ^ J, (17) 
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where the same symmetry relations in Eqs. (|ll[-|12[) relate the particle-hole compo- 
nents of the self-energy. 

The matrix u(p,p') represents the effective impurity potential for scattering 
of quasiparticles from p — > p' on the Fermi surface. In the simplest formulation 
of the HSM we neglect spin-flip scattering by the aerogel, assume that the aerogel 
is isotropic on the coherence length scale and retain only the isotropic scattering 
potential. The neglect of spin-flip scattering is justified for 3 He-aerogel with a small 
concentration of 4 He added to displace the 3 He solid layers. The resulting scatter- 
ing potential is then modeled by a single s-wave matrix element, u(p, p') = uqI. 
The assumption of isotropic scattering is more questionable. However, anisotropic 
scattering and orientational correlations of the anisotropic scattering centers repre- 
senting the silica strands can be incorporated into the HSM if needed. 

The Landau molecular-field self-energy 

SFL(P)= l E FL + h(p).^J' ™ 

is a functional of the diagonal components of the quasiclassical propagator. The 
scalar and vector components represent internal fields generated by the quasiparticle 
interactions through the disturbance of the quasiparticle spectrum in response to 
pairing correlations or external fields. The scalar molecular field is given by 

s -(p) = TY / J^ L A°(p-p') g(p'; e n ) , (19) 

where A s (p,p') is the spin-independent quasiparticle-quasiparticle scattering am- 
plitude. Similarly, the vector component of the molecular field, 

h(p)=T^y^^(p.p')g(p'; £ „), (20) 

represents the internal exchange field acting on quasiparticles at the position p on 
the Fermi surface, and is determined by the quasiparticle-quasiparticle exchange 
scattering amplitude, A a (p, p'). These forward scattering amplitudes are related to 
the Fermi-liquid parameters, i^ , through the Legendre expansion]!!! 

Ais ' a) = E jhr, p ^ ■ p') ■ ( 21 ) 

V l + FJ '/2£+l 

In calculating the linear response to a uniform magnetic field only the exchange 
interactions, Fq and F%, contribute for isotropic impurity scattering in the Born 
and unitarity scattering limits. 

Finally, we have the Zccman coupling to an external magnetic field, 

v mt = h cxt • ± , with h cxt = -Z a M H , (22) 

S = |(1 + T3)er + i(l — T3)<T tr , and [i — jH/2 is the nuclear magnetic moment of 
the 3 He nucleus. The factor Z§ = 1/(1 + Ffi ) is the high-energy renormalization 
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of the Zeeman coupling of the quasiparticle moment to the external field. Spin- 
rotation invariance in normal liquid 3 He fixea_the renormalization factor in terms 
of the quasiparticle exchange interaction, Ffi S3 

The quasiclassical transport equation and constitutive equations for the self- 
energies are_s.upplemented by Eilenberger's normalization condition for the matrix 
propagator ,£3 

9 2 = -tt 2 1. (23) 

In addition to the overall normalization of the propagator, the normalization con- 
dition is useful in simplifying the linearized transport equation for the response to 
an external field. 

4.1. Homogeneous equilibrium in zero field 

For the Balian-Werthamer phase the order parameter is given by 

A(p)=AR-p, (24) 

where R is the relative spin-orbit rotation matrix that defines the B-phase.Ej The 
BW state is isotropic under joint spin and orbital rotations; as a result the mag- 
nitude of the order parameter, |A(p)| = A(T), is also isotropic. The magnetic 
susceptibility is determined by exchange interactions and nonmagnetic scattering; 
we can safely neglect nuclear dipolar interactions, and for the purpose of calculating 
the susceptibility we can set R = 1 . 

The homogeneous equilibrium solution for the propagator and self-energies 
of the bulk BW phase is straight-forward to generalize for isotropic scattering, 
and has been imzjestigated by Buchholz and Zwicknagl for the BW model for a 
superconductor .Ej The molecular field corrections vanish in zero field, so we are left 
with the mean-field order parameter and impurity self energy in the homogeneous 
transport equation, 

ie n r 3 - E imPi0 - A , flo = , (25) 

where Qq is the equilibrium propagator in zero field. The solution for go for an 
isotropic BW phase with isotropic scattering has the same matrix form as that for 
the pure BW phase, 

ii n T 3 - A (p;e„) 
So = -71" ; — , (26) 

+ |Anl 2 



where the renormalized Matsubara frequencies and order parameter are related to 

zero-field impurity self energy by (we drop the subscript '0' unless it is explicitly 
needed) 

ii n = ie n - s 3 (£„) , (27) 

A(p;e„) = A(p) + A imp (e n ), (28) 
with the zero-field impurity self energy given by, 

£ imp = s x (e„) 1 + S 3 (e„) f 3 + A imp (e n ) . (29) 
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However, the off-diagonal impurity self energy, A imp (e„), vanishes for pure s-wave 
scattering because the the driving term in the equation for A(p;e„) is the mean 
field order parameter A(p), which vanishes when averaged over the Fermi surface,EJ 

(A(p))=|^A(p)=0. (30) 

Thus, A = A(p), and the t-matrix is isotropic and reduces to diagonal form, 

t = hi + hh, (3i) 



_ y/gj - a) el + A 2 
1 " 7rN f e^ + A^l-a)' [d) 



_ a ie n ^/e 2 + A 2 
3 ~ nN f el + A2(l - a) ' {AA) 

where a is the normalized scattering cross-section and is related to the scattering 
potential, uq, and scattering phase shift, <5o, by 

g-ato'ft- fryV (34) 

1 + {-KNfUoY 

The unitarity limit corresponds to strong scattering, u$ — ► oo, in which case a — > 1 
(5q — > 7r/2). Born scattering is the limit of small cross section, a < 1, (<5 < 1). 
The other key parameter needed to relate the HSM parameters to the properties 
characterizing the aerogel is the mean free path, £. In the HSM the mean-free path 
is determined by the Fermi velocity and the mean time for a quasiparticle to scatter 
off an impurity, I — v/t. The scattering time is determined by the imaginary part 
of the normal-state impurity self energy, 

„ h n, m „ 

Thus, the isotropic HSM is defined by two parameters which we take to be the mean 
free path, £, and identify with the geometric mean free path for a specific aerogel; 
and the dimensionless cross-section of the scattering centers, a, which can vary from 
(7 = (Born limit) to a — 1 (unitarity scattering limit). Since the geometric cross- 
section of a typical aerogel strand is large compared to the Fermi wavelength the 
unitarity limit is more appropriate for modeling the effects of scattering by aerogel. 
Finally, the renormalized Matsubara frequencies are given by the solution to 
the equation, 

£« - £n + A N - 2 , A2 /-. TT • (OO) 

el + A 2 {l-a) 

The equilibrium equations in zero field are closed by the self-consistency equation 
for the order parameter, 

^-^'ff^np-p')-^, 07) 
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which reduces to the BCS gap equation, 

A = V 1 ttT J2 7^ ^ - (38) 

^ y/el + A 2 

The linearized gap equation for pure 3 He-B [e n — > £ n ) determines the bulk transition 
temperature, 



k„|<£c 
l/tq = 7rT c0 J2 r^ = K(T c0 )~Hl.l3e c /T c0 ). (39) 



Using Eq. (k59J) we can remove the pairing interaction and cutoff in favor of the bulk 
transition temperature, T c q. Thus, the gap equation for 3 He-B-aerogel becomes, 

ln(T/T c0 ) =,TV [ -— L== - ±A . (40) 






V% 



The linearized gap equation reduces to the Abrikosov-Gackov formulae^ for the 
suppression of T c by scattering in the HSM for the aerogel,E3 



T N e 



ln(T c /T c0 ) = 2ttT c ^ ( 7—^ - — ) . (41) 



Density of States 

In superconductors with an order parameter that breaks orbital rotation sym- 
metry, such as d-wave pairing in the cuprates or f-wave pairing in UPt3, scattering 
from an impurity, particularly in the strong scattering limit, or from a surface, leads 
to the formation of quasiparticle states at the Fermi level that are bound to the 
impurity, or surface, within a distance of order the coherence length.EjEj In the 
case of a random distribution of scattering centers the impurity iupund states at the 
Fermi level broaden into an impurity band near the Fermi levelHa The same is true 
for the p-wave phases of 3 He in aerogel. 

The equilibrium solutions for the quasiclassical propagator also determine the 
quasiparticle excitation spectrum below T c . Continuation of the Matsubara Green's 
functions to real energies determines the retarded equilibrium Greens functions, 

g R (p;e)=Q{p;ie n ^e), (42) 

which provide spectral information for the low-energy quasiparticle states in the 
presence of pairing correlations and disorder. In particular the angle-resolved local 
density of states, averaged over spin projections, is given by 

N(p;e) = --Irag R (p;£), (43) 

whereas the spin-polarization spectral density is given by 

P(p; £ ) = --Im g ii (p; £ ). (44) 
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Fig. 3. The density of states for the BW state in homogeneous aerogel with a 
mean-free path of £ — 1800 A and cross-section of a = 0.99. The DOS spectra 
correspond to pressures of P — 5, 10, 15, 20, 25, 30 bar, and were calculated 
at the same reduced temperature for each pressure, T — 0.2T C . 



Figure shows the self-consistently determined density of states (DOS) for 
3 He-B-aerogel calculated for £ = 1800 A and a = 0.99 as a function of energy for 
pressures spanning the range P = 5 — 30 bar. The pairbreaking parameter 



£o 



2tttT, 



,i) 



(45) 



decreases from a = 0.23 — ► 0.09 over the same pressure range. Note that even at 
the highest pressure, where the pairbreaking effect is weakest, there is a substantial 
density of quasiparticle states at the Fermi level from the impurity-induced Andreev 
bound states. At lower pressures the B-phase becomes completely gapless. These 
quasiparticle states lead to an increase in the spin susceptibility. Their polariz- 
ability is described by the polarization spectral density. In particular, the change 
in magnetization can be expressed in terms of P(p;e) by analytic continuation of 
Eq.(||) to real energies, 



m = 2 



+ 00 



def(e)(P(p;e)), 



(46) 



where f(e) is the Fermi function. 

4.2. Linear Response to the Nuclear Zeeman coupling 

The linear response of 3 He-aerogel to the nuclear Zeeman field is obtained by 
linearizing the propagator in the external field, 

8 = io+|i + -». (47) 

and similarly for the order parameter and self-energies; Qq is the zero-field prop- 
agator given in Eq. (|2q). The first-order correction, $ji, satisfies the linearized 
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transport equation, 

ie n r 3 - s imp , - A , §1 = w cxt + Ex + Ai , go , (48) 

and the corresponding linearized normalization conditions, 

0o0i+0i0o = O, with qI = -tt 2 1. (49) 

The operators on the left-side of Eq. (ji^) can be expressed in terms of jjo> 

ie n f 3 - s imp . - A = AT -1 jj - n top ti 1 , (50) 



with N = -ir/yje 2 n + Ag. Thus, Eq. @ becomes, 



[So , 0ij 



N 



Ei+Ai, 0o 



(51) 



Combining Eq. ( plj ) with the normalization conditions in Eq. (fl9j) gives the linear 
correction to the matrix propagator in terms of the external field, first-order self- 
energies and zero-field propagator, 



27r V / ir 



A^ 



{**(< 



Ei+Ai 



) + 0o (i 



•Si 



Ai)0o} 



(52) 



The detailed solution for §i is most conveniently represented in terms of components 
of §i which are symmetrized between the particle and hole components, i.e. for any 
particular component of jji define X ± = X ± X. The symmetry relations (Eqs. 
[11-12 ) imply that X^ will have even or odd symmetry with respect to (p, e n ) — * 
(-p, -£„) or (p,e n ) -> (-p,e n ). 

The self-consistent solution for gi requires the first-order corrections for the 
internal field, the order parameter, and self-energies. We can express the correction 
to the impurity self-energy in terms of the first-order correction to the t-matrix, 



*i 



Nfi (Qi}i Q . 



(53) 



where (...) denotes the average over the Fermi surface. Equations (J31|-p3[) and 

( p2| ) for £ and 0i determine the first-order corrections to the impurity self-energy, 
which have the general structure given in Eq. (|17j). In terms of the symmetrized 
components, the linear corrections to the impurity self energy have the form, 



imp,l 


= 


"i 


*w 


= 


rii 


h+ . 

imp : l 


= 


"i 


h- p . 


= 


«i 


*£», 


= 


"i 


A± , 


= 


n. 



> P N f {(tt + 4)(gt) + 2tit 3 (gi)} , 

^/{(tf+ti) <sr> + 2*1*3 <^>} , 
v Nf {(* i +*i)<gi-) + 2t 1 t 3 <gr)} , 
»Nt {(*i+*i)<gr>- 2 *i*3 (si + )} , 

ip N f (tl-t 2 3 ){f±). 



(54) 

(55) 
(56) 
(57) 
(58) 
(59) 



The key corrections are the internal fields, h^ u and the impurity correction to 
the order parameter, A im ± . The scalar corrections, Hr^ i; turn out to vanish. 
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The external field also induces an internal exchange field given by Eq. (|2C|). 
Thus, the total effective field is 

n+(p; £rl ) = h+ pl (e n ) + h+ a (p) + 2h cxt , (60) 

h~(p;e„) = h i - pl (e„) + h FL1 (p), (61) 

with the exchange terms given by 

h F ± L , 1 (p) = /^f i ^ a (p-p')T^g 1 =(p'; e „). (62) 

In general there is also a linear correction to the spin-triplet order parameter, 

Af(P) = J^V t (p-p>)TY / iHP'^n), (63) 

which may be combined with the off-diagonal impurity correction, A 1 (p; e„) — 
^i (p) + Aim i( e «)- The full solution for gi can be expressed in terms of these 
effective fields: 

(64) 

(65) 
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= u, 
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ir (A A + N l 


9\ 


( e -2 +A 2 )3 / 2 ?C « [ A A lJ ' 


gi + = 


n A(A-h+) 


gi _ = 


-- ( e -2 +A 2 ) 3/ 2 { Ax(Axh ) + £ „AxA 1 } 


ft - 


= (£ -2 +A 2)3/2( -«) A h+ ' 


fr = 


= o, 


f! + = 


- *** n M*l z, + ) 




V4 + A2 ( £ ~2 +A 2 ) 3/2 y 1^- 


fl" " 


tAj 7T _ 

■ J frAxhlA 




V4 + A 2 ( £ l + A 2 ) 3/2i "" 



(66) 

(67) 
(68) 
(69) 
(70) 



We can further simplify the results by noting that A 1 obeys homogeneous, linear 
self-consistency equations which allow only the solution A. 1 = 0. Thus, g± and fj*" 
also vanish. 

The magnetization of 3 Hc-B in aerogel can be calculated from the Fermi surface 
average of Eq. ( pq ) and the corresponding total field. Thus, the susceptibility (Eq. 
(@)) is given by 

Note that only the even (+) part of g contributes to the magnetization. We can 
calculate the change in the magnetization, 

m = T^ fct) > (73) 
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from the solutions for gi* in Eqs. (|6q-p7|), and the self-consistency equations for 

h* - Eqs. (56 54|62|) - and A x given in Eq. ( |59| ) and Eq. (J7l|). Inserting the result 
for gj from Eq. (pq) gives, 



m = T^^+( £n )(p(p-h+)), 



(74) 



where 



K+(e n ) 



irA 2 



p2 
L°r- 



v2l3/2 



(75) 



The Fermi surface average of the external field and the exchange contributions to 
the total field can be evaluated, or expressed in terms of m, 



(p (p • h+)) = -h cxt + - A a r. m + — AS m + - h+ 



(76) 



We then obtain an equation for m that depends on the impurity contribution to 
the internal field, 



1 - l -A« ya - 1a$ yf) m = ~ 2/| h cxt + It J2 K + (s n ) h+ 



(77) 



where 



r]T /! + (£„) 



(78) 



is related to the Yosida function (Y = 1 — y 3 ) evaluated with renormalized Matsub- 
ara frequencies. To complete the calculation we need to calculate the sum over the 
impurity correction to the internal field, h+ . In general, for scattering that is nei- 
ther in the Born or unitarity limits, the self-consistent solution for the impurity field 
entering Eq. ([T7j) involves both the even- and odd components of the quasiparticle 
polarization, g 1 (see Eq. (poj)). However, for the Born and unitarity scattering 
limits the coupling to g^ drops out, and an analytic result can be obtained for 
h+ p (e„). 

In the unitarity scattering limit, a — ► 1, in which case t\ — ► and 



h 



Thus, the equation for hr^ (e„): 



1 v/fl+A^ 



irNf 



is, 1 



(79) 



1 

27TT 



K+(e n ) <p( P -h+)> 



(80) 



depends on the same Fermi surface average that determines m in Eq. (74), in which 
case we can solve for h+ , and then TJ2 e ^ + ( £ n)h^ n (£ n ), in terms of h oxt and 
m: 

I Z»'+ To / 1 o \ 

(81) 



2ttt 



(i^) + ek K+ 



3 



3 ° 15 2 
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The result for the change in magnetization reduces to 



m 



§ n h = 



l -\ A lvi-h A, iyi_ 

where the impurity-renormalized response function is given by 



(82) 



- r £^ 



1 



[el + A2]3/2 ] 1+JL (* 



(unitarity) . 



(83) 



Thus, the susceptibility of the BW phase of 3 He in aerogel in the HSM is obtained 
from Eqs. ©, (Frf) and ©, 



XB 

XN 



= 1- 



1 



\h 



i + ^oV \i-Uj% + &4)v 



(84) 



This result can be expressed in terms of the conventional Landau parameters, Ftf 



and F 2 °, 



Xb/x 



(1 + F$) 



Y(: 



;F 2 a ) 



N 



(85) 



1 + F Q »(§ + If) + 1^(1 + §f ) + l F a F af ' 

where Y — 1 — y 3 is the impurity renormalized Yosida function given in Eq. (M). 

In the Born limit, a — > 0, we can carry through essentially the same analysis 
as in the unitarity case to calculate the impurity contribution to the local field and 
obtain, 



1 



K* 



2lTT V 1 - 7T- K^ 

DTTT 



3 



3 u 15 2 



(86) 



For the Born limit the susceptibility is again given by (p4|), but with 

A 2 i I 



H 



^£71 



K + A 2 ] 



213/2 



1 



A 2 



(Born) , 



(87) 



6t [e=+A 



ji-s 2 



which yields Eq. (g). 



4.3. Ginzburg-Landau Limit 



For temperatures just below T c the mean-field spin susceptibility for non-ESP 
states decreases linearly with temperature. The temperature interval over which the 
susceptibility is linear defines the region of applicability of the Ginzburg-Landau 
(GL) theory. The general result for the susceptibility of the B phase in the GL 
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region is expressed in terms of the GL material parameter for the change in the 
Zeeman energy, i.e. Af2 Zcoman = g z A 2 H 2 jEi which implies 

Xb - Xn = -2 g z A(T) 2 ~ -g z T c (T c - T) . (88) 

We can obtain the GL limit from the general result for \b and identify the GL ma- 
terial parameter g z for the HSM of 3 He-aerogel. For T — ► T c , the spin susceptibility 
reduces to 

XN 3 \1 + F$J **' V ' 

with 

Thus, the spin susceptibility reduces to 

£-Ki**) (*£)'*«■ 

where 

S m (x) . £ (-^) , m>l, (92) 

n>0 2 ^ 

and x = 1/4tttT c is the pairbreaking parameter. The order parameter in the GL 
limit can also be expressed in terms of T c and the sums defined in Eq. (^2|) , 



1 - xS 2 (x) 



Y V^3(^) + (2CT-l)x5 4 (a;)y cV 7 v ' 

The result for the GL material parameter that determines the change in spin 
susceptibility is then identified as 

g z = ^ -^_ (94) 

y (1 + F a ) 2 127r 2 T c 2 ' K ' 

For a; = we recover the known weak-coupling result for g z for pure 3 Hc-B, 

_ _ ^/M 2 7C(3) 

ffz (1+^)2 12^^- ^ 

5. Susceptibility 

In Figure |J we show the impurity-renormalized Yosida function, Y at P = 
10 bar in the unitarity limit of the HSM for mean-free paths of £ = 1700 A, £ = 
5000 A and £ — > oo. The increase in F(T = 0) reflects the density of fermionic 
states below the gap that can contribute to the low-field magnetization. We also 
show the result for Y in the Born scattering limit at £ = 1700 A. The smaller cross- 
section is a less effective pairbreaker, which is reflected in the stronger condensate 
suppression of the spin-polarization response function. 



E-i 
l>- 




T/T 

' C.a 

Fig. 4. The renormalized Yosida function for the spin response function for 
the BW state vs. T/T ca for the unitarity limit (solid lines) of the HSM for 
P = 10 bar and I = 1700 A, t = 5000 A and £ = oo (bulk). The result for the 
Born limit (dashed line) is shown for t — 1700 A. 
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Fig. 5. The susceptibility for the BW state vs. T/T ca scaled to \n in 
the unitarity limit (solid lines) of the HSM for P = 10 bar and t = 1700 A, 
£ = 5000 A and £ = oo (bulk). The GL slope for I = 1700 A is plotted 
to indicate the limited range of validity of the GL result. Also shown for 
£ = 1700 A is the effect of "switching off" the impurity correction to the 
local field in the unitarity limit (dotted curve). The Born limit for \b is 
shown as the dashed curve for I — 1700 A. The material parameters used for 
P = lObar are: v f = 45.66 m/s, T c0 = 1.83 mK, F a = -0.744, F 2 a = -0.466. 



The magnetic susceptibility for 3 He-B in aerogel is shown in Fig. |5j for the same 
pressure, scattering cross-sections and mean-free paths. The Ginzburg-Landau limit 
for unitarity scattering and £ = 1700 A is also shown. The Fermi- liquid parameters 
are taken from the tables in the review article by Halperin and Varoquaux.o The 
effect of scattering by the aerogel is substantial; T c is suppressed to T c rs 0.5T c o 
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at P = 10, bar while the susceptibility is increased by nearly a factor of 2 at low 
temperatures. Another feature to note is that the Ginzburg-Landau limit for xb is 
valid only in a limited range of temperatures below T c ; thus, linear fits to limited 
data for x(T) may lead to misleading information on the material parameters when 
compared with the GL result. 
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Fig. 6. The density of states for the BW state vs. energy, e, for both 
the unitarity and Born scattering limits in the HSM for P — 10 bar and 
t = 1800 A. 







e/A 

Fig. 7. The spin polarization spectral function for the BW state vs. energy, 
e, for both the unitarity and Born scattering limits in the HSM for P — 10 bar 
and £ = 1800 A. 



Also shown in Fig. g for I = 1700 A is a comparison between the unitarity and Born 
limits for the low temperature susceptibility. Unitary scattering is more effective 
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at pairbreaking than Born scattering. In the same figure we show the effect of the 
impurity vertex correction to the internal field (Eq. |80| )) on the spin susceptibility 
at I = 1700 A in the unitarity limit; the impurity field leads to a relatively small 
correction compared to the overall increase in the susceptibility that results from 
the polarizability of the quasiparticle states. The large density of states near zero 
energy (see Fig. 0), characteristic of the unitarity limit, contributes relatively little 
to the polarization, as is evident from Fig. 0. As a result the quasiparticles states at 
finite energy below and near the gap edge, which are present for both unitarity and 
Born scattering are mainly responsible for the increase in the spin susceptibility. 
The integrated polarization from the gapless states below the Fermi level is larger 
in the unitarity limit reflecting the stronger pairbreaking effect for scattering from 
larger cross-section impurities. 

5.1. Comparison with Existing Data for the Susceptibility 



Experimental data £pr the nuclear spin susceptibility of 3 .He in aerogel is re- 
ported by Sprague, et al.cl for P = 18.7 bar and by Barker, et alB for P = 32 bar. In 
both experiments the susceptibility was reported to decrease below T ca when several 
monolayers of 4 He were added to suppress the solid-layer Curie susceptibility. 



p=18.7bar 



ftJ^t*V^ ^ 




T(mK) 

Fig. 8. Magnetic susceptibility of ^He in aerogel at P — 18.7 bar from 
Sprague, et al. The bulk susceptibility data (open circles) is also shown 
for the same pressure. The Landau parameter Fq — —0.757 is taken from 
the normal state susceptibility, while F£ is used to fit the bulk data to the 
theoretical result. The fit parameter is: F$ = 2.2. The theoretical results 
for 3 He-aerogel are shown for the same pressure and Landau parameters for 
mean free paths of £ = 1700 , 1800 , and 2000 A, and a = 1. 



In Fig. g we compare the susceptibility data at P = 18.7 bar with the theo- 
retical result calculated for the same pressure in the unitarity scattering limit. The 
comparison is only partly satisfactory. The pure He-B susceptibility reported for 
the same pressure by Sprague, et al.El does not agree precisely with the theoretical 
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result using the_Fermi-liquid parameters based on the susceptibility measurements 
of Hoyt, et al.EJ However, if we fix Ftf by the normal state susceptibility, use the 
Landau parameter F£ to fit the bulk susceptibility data shown in Fig. g, calculate 
the susceptibility for 3 He-aerogel with these Landau parameters, then the compar- 
ison with the data for 3 He-aerogel is qualitatively correct for reasonable values of 
the mean free path. 

The comparison between the HSM for 3 He-B in aerogel and the data of Barker 
et al. is slightly improved. The theoretical calculation shown in Fig. was ob- 
tained with the bulk Fermi-liquid data from Halperin and Varoquaux's tables, and 
a mean-free path of i = 1800 A for unitarity scattering. The data of Barker, et al. 
agrees semi-quantitatively with the theoretical result, but the data span a limited 
temperature range and shows an unphysical upturn at lower temperatures, which 
is also visible in the data of Sprague, et al. even for the bulk measurements. 



15 



a 

so 

a 
a 



10 



p 


i i i i i i 
= 32bar 


l i i i i l i i 


, 


1 = 


= 180nm 




- 


a 


-+1 






- 




■TT 


jF** 














%/ 




- 


i 


!.,,.!,, 


, , i , , 



0.5 



1.5 



T(mK) 



Fig. 9. Magnetic susceptibility of 3 He in aerogel at P = 32 bar from Barker, 
et al. The theoretical result calculated for the same pressure, a mean free 
path of £ = 1800 A for 9 = 1. 



6. Conclusion 



Experiments based on NMR and acoustic attenuation indicate that the equilib- 
rium phase of superfluid 3 He-aerogel at all pressures is a non-ESP phase, generally 
assumed to be the Balian-Werthamer state modified by scattering off the aerogel 
structure. For 3 He-aerogel impurity scattering leads to substantial changes in the 
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susceptibility of the BW phase. The available data is in semi-quantitative agree- 
ment with theoretical predictions. Additional measurements could prove important 
for a more definitive identification of the equilibrium order parameter, as well as for 
refining the theoretical model to describe the effects of disorder and scattering on 
the properties of superfluid 3 He. 
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